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We report on the geometry and mechanics of knotted stiff strings. We discuss both closed and 
open knots. Our two main results are: (i) Their equilibrium energy as well as the equilibrium tension 
for open knots depend on the type of knot as the square of the bridge number; (ii) Braid localization 
is found to be a general feature of stiff strings entanglements, while angles and knot localization 
are forbidden. Moreover, we identify a family of knots for which the equilibrium shape is a circular 
braid. Two other equilibrium shapes are found from Monte Carlo simulations. These three shapes 
are confirmed by rudimentary experiments. Our approach is also extended to the problem of the 
minimization of the length of a knotted string with a maximum allowed curvature. 



I. INTRODUCTION 

We report on the properties of stiff knots, i.e. knotted 
strings whose shape is dictated by the bending curvature 
energy and contact interactions only. Stiff knots, such 
as loose knots with nylon or metal strings, are ordinary 
objects in everyday life (Fig.l). An upsurge of interest 
in stiff knots recently came from studies which pertain to 
microscopic objects, such as those encountered in biology 
and nano-technologies. Some examples are: knots with 
actin filaments^, nanotubes^, nanotube fibers^, and silica 
wires^ (see Fig. 2J. These studies point out the relevance 
of stiff knots for the experimental determination of the 
bending rigidity 1 , for knot induced polymer and filament 
break-upi^, and for nano-manipulation 2 . 

Moreover, many recent experimental and theoretical 
studies were devoted to knots in polymers, with an em- 
phasis on knotted DNA 6 . They show that flexible poly- 
mers are subject to knot localisation leading to the forma- 
tion of small prime knots. Knot localization may result 
from entropic effects^, or long range interactionsii*i£. If 
the localized knots are small enough, their thermal fluc- 
tuations become negligible and they might be described 
by the stiff knot regime. 

Finally, stiff knots may be considered as elementary en- 
tanglements which capture qualitatively some of the fea- 
tures of more complex entanglements. Hence, stiff knots 
may also provide insights for the curvature energy dom- 
inated behavior of tightly entangled semi-flexible poly- 
mers such as acting, and other fibrous materials^. 

Here, we aim to establish the basic mechanical and ge- 
ometrical properties of stiff knots. First, the mechanical 
properties of the knots are found to exhibit a surprisingly 
simple dependence on the knot type via a quantity called 
the bridge number n (to be defined below). We show that 
the minimum knot energy, as well as the minimum equi- 
librium tension at the open ends of an open knot, both 
increase with n as n 2 . Secondly, we identify a striking 
and general feature of the geometry of stiff strings entan- 
glements, which we call braid localization. We analyze 
the geometries of the simplest knots which indeed exhibit 
braid localization. 

We will begin with a study of closed knots. In section 




FIG. 1: Stiff knots. Upper panel: A closed knot, made with 
a plastic tube of width 1 cm and length 80 cm (see section 
IIX CI for details). Lower panel: Open trefoil knot with a silica 
nanowire, from Ref^. The width of the picture is 20 /im. 



II, we define the curvature energy of a filament. In section 

III, we discuss the role of interactions on the equilibrium 
shape of a filament, and we place our problem within 
the historical perspective of the studies on the so-called 
Bernoulli-Euler elastica. A lower bound for the energy 
of the equilibrium configuration is given in section IV. 
In section V, we analyze the limit of thin strings which 
leads to braid localization. In section VI, the results of 
section V are used in order to obtain upper bounds for 
the equilibrium energy. Two results follow from the exis- 
tence of these upper bounds: (i) the global minimum of 
the energy scales as n 2 ; (ii) the minimization problem is 
solved exactly for a special family of knots. The shape 
and energy of some simple knots are obtained from Monte 
Carlo Simulations in section VII. We find three geome- 
tries which correspond to the simplest knots. In section 
VIII, we translate the main results of our analysis to the 
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case of open knots. In section IX, we discuss several ad- 
ditional points: (i) the equivalence between the limits of 
thin strings and long strings; (ii) the curvature energy 
of thick knots; (iii) rudimentary experiments; (iv) the re- 
stricted curvature model, which makes link with the work 
of Buck and Rawdon 27 . Finally, we conclude in section 
X. 



II. MODEL 

We shall first focus on closed knots, i.e. single knotted 
strings without free end, as in the upper panel of Fig^ 
A given configuration of a knot in 3D space is described 
by a position vector r(s), where s is the arclength. Since 
the knot is closed, r(s) is periodic in s, and its period is 
the length 




of the knot. In the following, the absence of bounds in 
the integrals indicates integration over the whole knot. 
We define the usual curvature energy of an inextensible 
string 1 ^^ as: 

E = ^JdsK 2 : (2) 

where k > is the curvature, and C is the bending rigid- 
ity. Such a modelling is valid in the limit of small defor- 
mations, defined by the limit of small curvature^: 

WK«1 (3) 

where w is the diameter of the section of the filament. 

The question we address is the following: for a given 
knot K of length C, what are the equilibrium shape and 
energy? We call equilibrium energy the global mini- 
mum energy, as opposed to that corresponding to pos- 
sible other local minima, which will be referred to as 
metastable states. The equilibrium energy is denoted as 
£* . We shall fix C by means of a Lagrange multiplier /i. 
Mechanical equilibrium is then obtained from the mini- 
mization of: 

T = £ + fiC=j I dsn 2 +fi I ds. (4) 



III. ON THE ROLE OF INTERACTIONS 

Let us first consider this minimization problem in the 
absence of interactions between the different parts of the 
strings. In this case, the curve can freely cross itself. 
We obtain the so-called elastica, initially proposed by D. 
Bernoulli. In order to investigate further this problem, 
we consider the variation of T induced by the variation 



Sr(s) of the position of the string. Since the filament is 
closed, there is no boundary terms and: 

5T = J ds Sr.d s A (5) 

where 

A= ^K 2 -/ijt + Ca s Kn + CKrb. (6) 

Throughout the paper, we denote the derivatives with 
respect to s as 9 S . At equilibrium, the energy variation 
vanishes: 5T = 0, leading to the nonlinear differential 
system d s A = 0, i.e. A is constant, and^ 

|« 2 - M = A.t, (7) 
-Cut = A.b. (8) 

The constant vector A represents the internal forces in 
the string^, r is the torsion, t = d s r is the tangent vector 
of the curve, and (t, n, b) is the usual Frenet frame. We 
have not written down the projection of A on n because 
it is redundant. Indeed, it and can be obtained from a 
derivation of (Q with respect to s. 

But 5J- does not vanish only at the global minimum 
(i.e. at equilibrium), and a number of other spurious 
solutions are found? 8 . Planar solutions of (|7I8() were an- 
alyzed by Eulenji^. The 3D solutions are listed in RefE 
Comparing their energies, one finds that the circle is the 
closed solution with the lowest energy. It is also the only 
closed solution which is stablsi^. Thus, knots cannot be 
stable solutions of I7I8|I . Knots can nevertheless be stabi- 
lized in the presence of interactions between the different 
parts of the string. 

In the following, we use a hard core repulsion, and 
the string is modelled by means of a non-self-intersecting 
tube of diameter w. Such a hard core repulsion implies 
that the distance between different parts of the string is 
larger than w, and the radius of curvature 1/k is larger 
than w/2 (see e.g. Ref. IT^) . In the previous paragraph, 
we concluded that knots could not be stabilized in the 
absence of interactions. In the case of hard core repul- 
sion, such a statement means that contact points must 
be present. 

But the non-crossing condition at contact points in- 
volves interactions between distant parts of the string 
(i.e. parts with different values of s). These interactions 
can be included in the energy J2J) by means of an addi- 
tional term which includes a nonlocal interaction poten- 
tial. In the variational formulation, such a term leads to 
additional nonlinear integro-differential contributions to 
(|7I8|) . the consequences of which are difficult to analyze. 
Despite the difficulty, some mathematical informations 
about the solutions - such as their existence- were ob- 
tained from this approach 16 . 

We here consider a different approach based on a com- 
bined analysis of the knot topology and of the geometry 
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(a) (b) 

FIG. 2: The trefoil knot can be represented in different ways 
in 3D space such as in (a) and (b). For each configuration, 
there is a given number of maxima along the x axis: 2 in (a), 
and 3 in (b). The minimum number of maxima among all 
configurations of the knot is the bridge number n. On the 
figure, we have drawn the trefoil knot, for which n = 2. 

in the limit of vanishing string width w — > 0. Our ap- 
proach does not systematically provide the equilibrium 
configuration and its energy, but it allows one to ob- 
tain important information about the equilibrium energy 
(such as lower and upper bounds), and about the geome- 
try of the equilibrium configuration (such as the absence 
of knot localization and the presence of braid localiza- 
tion) . 



IV. LOWER BOUND FOR THE ENERGY 

We start with a result found by J.W. Milnor 17 : for any 
given knot K, 

2irn<R, (9) 

where 

R = J dsn (10) 

is a dimensionless quantity called the total curvature^, 
and n is the bridge number of the knot K. To define n, 
let us consider a given direction in space x. As shown in 
Fig. [21 any given configuration r(s) of K in space has a 
well defined number of maxima along x. The minimum 
number of maxima among all configurations of K is n. 
Deforming the knot so as to place all maxima and all 
minima in two parallel planes of fixed abscissa along x, 
we obtain the linear braid configuration of figure^. The 
minimum number of loops at the top of the braid is n. For 
an unknotted closed string (usually called the unknot), 
n = 1, and for any other knot n > 2. For example, n = 2 
for most DNA knotsSl 

We define the normalized energy 



which depends on the knot shape, but is independent of 
the knot size and bending rigidity. Hence, e only depends 
on w/C and on the type of knot. Using the Schwarz 




FIG. 3: Knots in ordinary 3D space. Any knot (a) can 
be deformed to a braid and loop configuration (b) with n 
maxima, or to a closed braid configuration (c) with i strands 
(here n = 2 and i = 3). The lateral shrinking of (b) leads to 
the configuration (d) with the lowest total curvature R — 2nn. 
When the central braid of (b) is shrinked to a point, and the 
loops are solutions of Eqs. 171811 . we obtain the Point Braid and 
Loops configuration (e). In (f), the closed braid is laterally 
shrinked, and has a circular shape. 

inequality: 




and we find: 

2ttV < R 2 /2 < e. (13) 

This generalizes a result of Refill which reads in our 
notations: e > 87r 2 , and which can be obtained from 1)1 3|) 
with the additional knotting condition n > 2. 

Note that @ is sharp: for any knot, it becomes an 
equality for the laterally shrinked configuration of Fig. 2d. 
Nevertheless, (|13[1 is not necessarily sharp: it only pro- 
vides a lower bound for the energy. Since it is a lower 
bound for all configurations, it is also a lower bound for 
the equilibrium energy. This lower bound depends on the 
knot type via n and but does not depend on the filament 
width w. 



V. ASYMPTOTICS OF THIN STRINGS AND 
BRAID LOCALIZATION 

A. Variation of the energy with w 

The energy £^ of the equilibrium configuration of a 
knot K is a monotonically increasing function of w. The 
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proof of this statement is reported in Appendix A. We 
may therefore write: 

d w £* w > (14) 

Hence, if we consider the equilibrium shapes of a given 
knot K as a function of w, the lowest value of £^ will 
be reached in the limit w — > 0. On the opposite, the 
highest value of will be reached by the tight kno^S, 
also called ideal knot^S, which is the configuration with 
the highest possible value Wid of w authorized by the hard 
core repulsion. This can be summarized as: 

r<s* w <r md . (is) 

where £* = lim w ^o ££,. 

Here, we will not analyze the full dependence of £^ 
with respect to w, but we rather focus of the limit w — * 0. 
The reason of our focus on this limit is self-consistency. 
Indeed, the expression of the energy £ , defined in Eq.(J5J, 
is valid in the limit |(3} . Integrating 1(3} along the knot, 
we find wk <C C. Then, using J[J}, we obtain 



Since we consider a given knot (i.e. fixed n) with a fixed 
length (i.e. fixed £), the limit w —> is required for the 
energy £ w of a knot to be well-defined. 

B. Configurations in the limit w — > 

We now analyze the knot configurations in the limit 
w — > 0. They are obtained from a procedure in 2 steps, 
(i) Identification a possible structure {S} of the knot 
when w — > 0. (ii) Variational approach on the structure 
{5}. 

_Z. Structure 

In step (i) , we take the limit u> — > which may lead to 
"wild knots" , exhibiting knot accumulations or singulari- 
ties. We shall determine which types of accumulations or 
singularities are allowed, and which ones are forbidden. 

The 4 possible types of singularities or accumulations 
are given on Fig0] Angular points -as in FigHJi, and knot 
localization -as in FigQJ), are forbidden because they lead 
to an infinite energy. A rigorous proof of this statement 
is given in Appendix [Bj Here, we only provide an in- 
tuitive explanation: angles and knot localization may be 
obtained by decreasing the length Ci oc of a part of a curve 
to zero, keeping its shape fixed. An example of such a 
shrinkage is shown on Fig. where the size of the dashed 
box containing the angle decreases. Using Ea. (|ll|) . the 
energy of this part behaves as £i oc = Ce/Ci oc where e de- 
pends on the shape but not on the scale. Since Ci oc — > 0, 
one has £\ oc — > oo. 




FIG. 4: In the limit where the string width w vanishes, var- 
ious singularities or accumulations may occur. The thickness 
of the lines are used to represent various filament widths w. 
(a) angle. The dashed line defines a box whose size is shrinked 
to zero to obtain localization. The filament length in the box 
is Cioc- see text; (b) knot localization; (c) multiple point; (d) 
braid localization, leading to a line-braid, or a point braid. 
Angles and knot localization are forbidden because they lead 
to divergence of the energy £. 

The two types of accumulations which are allowed 
when w — > are multiple points -as in Fig. ^i, or braid 
localization -as in FigQJi, because they do not lead to a 
divergence of £ . 

Braid localization is defined as the lateral shrinking of 
a braid, with all strings in the braid tending to the same 
curve. We define a line-braid of multiplicity to > 1 as 
the result of the lateral shrinking of a braid of to strings. 
In the following, a single string will be considered as a 
line- braid with multiplicity m — 1. Since all strings in 
the braid tend to the same curve, the curvature energy 
of the line-braid reads mC/2 J ds k 2 . 

In some cases, we might then also shrink the length 
of the braid to zero to obtain a point-braid, as in Fig. 
0Jl. The order of the limits is important: firstly lateral 
shrinking, and secondly length shrinking, so that the cur- 
vature energy of the point-braid vanishes. Since angle lo- 
calization is forbidden, strings or line-braids must emerge 
tangentially along the same axis from a point-braid. 

By shrinking braids in a special way with a given knot, 
we obtain a structure {S} of line-braids connected by 
multiple points or point-braids. An example of such a 
structure is depicted on Fig[5J 

At this point, we shall put some emphasis on a central 
result of the present section: in the limit w — > 0, knot lo- 
calization is forbidden and braid localization is expected. 
We shall also stress on the fact that we have not proved 
that braid localization will occur, but we have shown that 
it can occur when w — > 0. 



2. Variational analysis of a structure 

In this section, we derive the conditions which must be 
satisfied at equilibrium from a variational approach. We 
consider a given structure {S}, with N line-braids. Its 
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m=l 




Point-Braid 



Line-braid 
m=2 



Multiple point 



FIG. 5: Shape of a knot in the limit w — > 0: an ensemble of 
line-braids connected by multiple points or point-braids. The 
drawing is 2D but we consider 3D configurations throughout 
the manuscript. 



energy reads: 



'{5} 



N 

x -s TUpC 



p=l 



dsn 



(17) 



where the index p in the integral means that the integra- 
tion is performed over the pth line-braid only, and m p is 
the multiplicity of the pth line-braid. As in sectional the 
total length is fixed by means of a Lagrange multiplier /i, 
and we have to minimize: 



{S} 



N 

E 

P =i 



m p C 



[im v I ds 



(18) 



The variation of 3~{s} resulting from a variation 8r(s) of 
the structure {S} reads: 



sjr {sy = E ( / ds Sr - d * A p 

P =i ^ J p 

+ [Sr.A p + Cm p Kn.d s 8r] p j (19) 

where [Y] p indicates the difference between Y at the end 
of the pth line braid, and Y at its beginning. Moreover, 



C 2 

2 K 



fj,\t + Cd s Kii + Cut b 



(20) 



At equilibrium, one has 5 J- = 0, which leads to d s A. p = 0, 
so that Ap is a constant vector. Projecting A p on t and 
b, we find: 



C 2 
— K 

2 



— Ap.t , 



-m p Cnr = A p .b. 



(21) 
(22) 



At equilibrium, line-braids therefore obey differential 
equations similar to that of single strings Q7I8[I . 

Let us now discuss the boundary conditions at the con- 
tact points (multiple points or point-braids) between the 



line- braids. We use the index B to list the points at 
which the line-braids are connected. At a point B, several 
strings related to different line-braids and point-braids 
may cross. Two types of constraints may force the strings 
which cross at B to be tangent to each other, (i) Since 
angle localization is forbidden, each string enters and ex- 
its at the point B with the same tangent vector, (ii) 
Each string is also tangent to some other strings because 
it belongs to a line-braid or a point-braid. The combina- 
tion of these two constraints forces a bunch of line-braids 
ending at B to be tangent to each other at the contact 
point. There might be several bunches of line-braids at 
the point B. These bunches can rotate freely from each 
other, but all braids inside the bunch are tangent to each 
other at the point B. Let qs be the index which lists the 
bunches at the point B. The boundary contribution of 
the variation l|19fl may then be re- written as: 

AT 

/J[— Sr.A p + Cm p Kn.d s 5r] p 

p=i 

= -E^- E A p 

B p^B 

+C^2^2d s Sr qB . m p Kn, (23) 

B qs V^IB 

where p — > B indicate that we perform the sum over all 
parts connected to B. Moreover, p S pb indicates that 
we perform the sum over all line-braids and point-braids 
belonging to the bunch qs- For dcfinitcncss, each braid 
should be oriented, so that each term in the sums over 
q — > B and p G qB contains a sign zbS. The quanti- 
ties Stb, and d s 8v qB respectively account for infinitesi- 
mal translation of the point B and rotation of the bunch 
qB around the point B. At equilibrium, one must have 
8F{S\ — f° r a U perturbations, so that from Eas. ijP9"|) 
and 



E a p = °- 

p->B 

m p n p n p = . 



(24) 
(25) 



p€q B 



The differential system (|21I22|) . together with the bound- 
ary conditions (|24I25|) is well posed and should be solved 
in order to determine the configurations of a given struc- 
ture {S} for which the variation of the energy vanishes, 
i.e. 6J-{s] = 0. In the following, we will denote the solu- 
tions which obey the equation 8J-^y = with the index 
°. For example, the resulting energy of a structure {S} 
will be written as £{gy 

As in the case of the Euler-Bcrnoulli clastica, which 
was discussed in section IIIII the class of configurations 
which obey the equation SJ^^gy = contains stable and 
unstable solutions. Our goal here is not to analyze this 
class of configurations in details. We will rather look for 
some specific configurations belonging to it, which will 
provide us with upper bounds for the equilibrium energy. 
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We shall now present two relations which apply to this 
class of configurations. The first one provides the general 
expression of the Lagrange multiplier M{s} f° r a given 
structure {S}: 



£ 



{S} 



C 



(26) 



The demonstration of this relation is reported in Ap- 
pendix |0 The relation l|41|) allows one to interpret fi° as 
the average curvature energy density (per unit length of 
string). 

Another relation, which applies to a restricted class of 
structures, will be very usefully in the following. Indeed, 
in some cases, the normalized energy of the structure 
can be expressed by means of the normalized energies its 
line-braids: 



N 
Lp=i 



-, 2 



Dl/2 



(27) 



and the length of the pth line-braid reads 



1/2 



(28) 



The situations where these formulas applies are: (i) all 
line-braids are arcs of circles, (ii) all line-braids are closed 
loops, and (iii) all line-braids have the same energy and 
the same length (this includes the case of line-braids with 
identical shapes). The derivation of this relation, as well 
as the general expression of e^ S y are given in Appendix 



VI. UPPER BOUNDS FOR THE ENERGY 

A. Bridge number and braid index 

Let us consider a first example of configuration belong- 
ing to the above-mentioned class. We have seen in section 
IIVI that any knot can be deformed in order to obtain a 
configuration similar to that of Figl^D with n maxima. 
In the limit w — > 0, this configuration can be deformed 
in such a way to obtain the Point Braid and Loops (PBL) 
configuration of Fig.2e, which defines the structure {S}. 
Each loop has the same boundary conditions: the curve 
begins and ends at the same point, the initial, and final 
tangent vectors being opposite. The detailed calculation 
of the minimum energy of one loop is reported in Ap- 



pendix El We find 6j° oop w 18.19. Since the structure 
{S} is an ensemble of identical loops, formulas I|27I28|) 
apply. Therefore, all loops have the same length and: 



o a 2 o 

e PBL — 4n e loop- 



(29) 



This expression provides a first upper bound for the equi- 
librium energy. 



A second upper bound is found using Alexander's 
theorem^!, which stipulates that any knot can be trans- 
formed into a closed braid, as shown on Fig[3J:. The 
smallest possible number of strings in the closed braid is 
the braid index and i > n. In the limit w — * 0, we 
may laterally shrink the closed braid with i strings, and 
we obtain a closed line-braid of multiplicity i. As men- 
tioned in section ITTT1 the equilibrium energy for a closed 
line is reached by the circular configuration, for which 
e, 

again obtained from 127|) 



o = 2ir 2 . The value of e° CB for a closed line- braid is once 



t CB 



I fcn 



2tt 



2 i 2 . 



(30) 



We see that for small i, i.e. when n < i < an, with a = 
{2ei oop ) 1 / 2 /n w 1.92, the CB configuration has a lower 
energy than the PBL configuration. On the opposite, for 
large i, i.e. when i > an, the PBL configuration has 
the lowest energy. Combining the lower bound ((13(1 and 
the upper bounds I(29I3U|) we find that the normalized 
equilibrium energy obeys: 



27rV < e* 



^ o 2 • r 2 2 

< 27T mm q n ; 



Using we finally have: 

022^/0*^-02 •r22 -2l^' 

L mm [a n ; % \ — 



(31) 



(32) 



These inequalities are a central statement of the present 
paper. Let us now present two results which directly 
follow from them. 



B. Scaling of the equilibrium energy 

First, the inequalities 1)32(1 allows us to reach a gen- 
eral conclusion: the equilibrium energy exhibits upper 
and lower bounds both proportional to n 2 C/C. We shall 
write this result as: 



e* 



(33) 



Although this relation does not provide the exact value 
of £*, it is a strong indication of the behavior of £* as a 
function of knot complexity. For example, we may con- 
clude that £* can be large only for knots having a large 
value of n. 



C. The n — i knot family 

The second consequence of (|32|l points out a specific 
class of knots. Indeed, when n = i, 1(3 2|l becomes an 
equality, and our problem is readily solved: we have 



£* 



2n 2 i 2 



C 
1 



(34) 



and the configuration is a circular line-braid of multiplic- 
ity i. The relation n = i defines a knot family which 
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contains the torus knots (knots obtained by wrapping a 
string on the surface of a torus without crossing). Us- 
ing available tabled we have analyzed prime knots of 
crossing number n c < 11. The crossing number n c is de- 
fined as the smallest number of crossings in the planar 
projections of a knot. We find that the n = i knot family 
contains w 20% of these knots. 



VII. MONTE CARLO SIMULATIONS 

In order to gain more insights about the configura- 
tion of stiff knots, we have performed Monte Carlo (MC) 
simulations with a closed chain of TV = 150 beads sep- 
arated by segments of fixed length -equal to 1. The 
length-preserving elementary motion of the chain is im- 
plemented via the rotation (with angle ±7r/100) of one 
bead around the axis which runs through its neighbors. 
We use the Metropolis algorithm, with the energy: 

N 

«d = Cj](l-u n .u B+ i) ) (35) 

n=l 

where u n is a unit vector along the nth segment. The 
closure of the chain imposes: 

N 

U n = , (36) 

n=l 

un+i = ui. (37) 

At low temperatures, the chain length C — TV is much 
smaller than the persistence length L p = C/ksT. Then, 
the curve is smooth, and Ed — > E. Non-crossing con- 
ditions are imposed with spheres of excluded volumes 
around each bead: we forbid beads to get closer than 
I/a/2. This leads to an excluded volume tube with a 
non-constant width w, varying between 1 and v2. We 
will analyze knots with n = 2 only, so that the condition 
(Q2J is verified: w/L ~ 10~ 2 < 1/2™ - 10 -1 . 

We use a simulated annealing method with a power law 
decrease of the temperature up to the low temperature 
regime. Repeated simulations with the same knot pro- 
vide us with the ground state and sometimes also with 
metastable states. 

A list of the knots studied in the simulations, to- 
gether with the obtained equilibrium configurations, is 
presented in FigEl A knot is usually denoted with its 
crossing number n c , and with an index which indicates 
its order in the standard list of knots for a given value of 
n c . For the trefoil knot -denoted as 3i, which is a torus 
knot with n — i = 2, the expected circular line-braid con- 
figuration is found, as shown on FigEl FigElalso shows 
two other configurations: the 4i (figure 8 knot), and 52 
knots lead to a configuration which will be denoted as the 
4- fold cage configuration in the following. The 52, 62, 63 
knots lead to a configuration which we call the 80 con- 
figuration (because it is composed of an 8 and a circle). 
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FIG. 6: Table of results for the MC simulations of closed 
stiff knots. Three different configurations are found: (i) the 
circular braid for the 3i knot; (ii) the 4-fold cage for the 4i 
and 52 knots (the sphere is a guide for the eye, indicating 
that this configuration is approximately wrapped around a 
sphere); (iii) the 80 for the 52, 62 and 63. 



The three different shapes are summarized on Fig[7| We 
observe that braid localization is present in all of them: 
line-braid for the n = i configuration, and point braids 
for the two other configurations. 

In Fig|S| we have shown the conjectured knot families 
which lead to the 4-fold cage and the 80. These con- 
jectured families are also reported on Fig|S] From knot 
tables—, the sum of the 3 families n = i, 4-fold cage, and 
80 represent most of simple knots: all prime knots with 
n c < 7, and 75% of those with n c < 9. 

We have analyzed in details the energy of the 80 con- 
figuration. It is calculated from the procedure of sec- 
tior JV Bl The structure of the 80 is composed of two 
parts: the 8, and the circle 0. The energy of the 8 curve 
is calculated in Appendix iDl eg ps 55.93. The is a cir- 
cle, and its energy is 6q = 27r 2 . Since the and the 8 are 
independent planar solutions, we may assume vanishing 
forces at their contact points, and analyze the structure 
of the 80 as a sum of two closed solutions (the and the 
8) . Using Eq. I|27|) , we may then evaluate the equilibrium 
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FIG. 7: (a) The 3i provides an example of line-braid of 
multiplicity 2. (b) The two contact zones in the 4-fold cage 
configuration are point-braids, indicated by white boxes, (c) 
An 80 configuration, exhibiting one point-braid and 3 double 
points, indicated by small circles. 




FIG. 8: Knot families leading to 4- fold cage and 80. The 
boxes account for twists with an arbitrary number of turns. 



energy of the 80 structure as: 



l 8<) 



- (4 1/2 



=1/2x2 
1 ) 



142.1. 



(38) 



As expected from section IV Al the energy in the MC 
simulations (given on FigEJ ; corresponding to a finite 
(but small) value of the width w, is slightly larger than 
the theoretical result corresponding to the limit w — ► 0. 

From the MC simulations, a normalized energy equal 
to 129.6 is obtained for the 4-fold cage. Since d w £^ > 0, 
this means that in the limit w — > 0, t\ c < 129.6. The 
value of e| c is probably quite close to its numerical upper 
bound 129.6, but we do not know its exact value. We 
only know from (|31|) that, since n = 2, one has t\ c > 
8ir 2 . Finally, the different configurations relevant for the 
simulations fall into the hierarchy: 



e Cs|i=2 < £4C < e 80 < e CBH=3 < CpBLl 



(39) 



As seen from FiglSJ some knots may exhibit more than 
one configuration. For example, both the 4-fold cage 
and the 80 configurations where observed with the 52 
knot. The hierarchy l|39|l shows that the 4- fold cage is 
the ground state, while the 80 is metastable. This result 
is also observed in the simulations: although the energies 
of both configurations are slightly larger than the asymp- 
totic value for w — > 0, their order is not affected by the 
finitcness of the width w. 



VIII. OPEN KNOTS 

Let us now consider an open knot, which exhibits two 
asymptotically straight open ends far from the knot as in 
Fig.l. We may ask the same question as for close knots: 




FIG. 9: Schematics of the open trefoil stiff knot. Open ends 
are along the z axis. From left to right, the filament width w 
decreases. 



what are the configuration and the energy at equilibrium? 
We shall here show that our results on closed knots can 
be extended to the case of open knots. 

For an open knot, the total length of the string di- 
verges. It is therefore more suitable to use the excess 
length C e , defined as the total length minus the length 
of a straight string: 



Le.= 



J ds(l - t ■ z) = J ds- J d: 



(40) 



where the integrals are performed along the whole knot, 
and z is the direction of the open ends far from the knot 
(We assume that both ends have the same direction). 
The excess length C e should replace £ in the expression 
of the energy T defined in Eq.Q. The variations are 
performed in the same way, and since the last term of 
(|40|l is constant, its variation vanishes, and the obtained 
differential system is the same as that for close knots. 

Since open ends are straight, we assume that k — > 0, 
and 8 s k — > far from the knot. Then, from (|20|l we have 
A.t — > fi° , A.n — > 0, and A.b — ► 0. Since A describes 
internal forces in the string, and t is the tangent vector, 
this means that ji° is the equilibrium tension far from 
the knot. 

For a given structure {S} of an open knot at equilib- 
rium, we show in Appendix [U] that the Lagrange multi- 
plier fi° is related to the curvature energy: 

u° - £ ° {s} -f° G UU 



(42) 



where the normalized energy of an open knot is defined 
as: 

o _ £ 

The general expression of e^fs} ^ s g lven m Appendix [O 
Note that relation l|41|l is also valid for closed knots, as 
shown in section IV B 21 Since the excess length C e is 
fixed, Eq. I|41l) also shows that the values of fi° are ordered 
in the same way as the energies. 

For open knots, the bridge number n\ and the braid 
index i\ are defined in the same way as for closed knots. 
If an open knot is obtained by opening a closed knot with 
bridge number n and braid index i, the bridge number 
and braid index of the open knots areS: 

ni = n - 1 , (43) 
h=i-l. (44) 
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As an example, since n = i = 2 for the trefoil knot, one 
has rix = i% = 1 for the open trefoil knot. 

The PBL configuration of Fig|3^ with open ends con- 
nected to the knot via a straight line tangent to the point 
braid, is used as a first upper bound. A second upper 
bound is obtained with the circular braid configuration 
of FigOF with open ends connected to the knot via a 
straight line tangent to the circular line-braid. Finally, 
(|32|l and i|41|) are also valid for open knots when C is re- 
placed by the excess length C e , so that the equilibrium 
tension /i* of open knots obeys 

< M* < 2TT 2 min[anl;i 2 1 }^ . (45) 

We therefore conclude that the equilibrium tension be- 
haves as nf, i.e. fx* ~ n\C/L 2 . Furthermore, when 
m = i\, the equilibrium shape of open knots is a cir- 
cle tangent to a point braid, and fi* = 2ir 2 i\C / C 2 . 

The simplest example of open knot is the open tre- 
foil knot (3i), for which n\ = i\ = 1. As depicted in 
FigHH our theory predicts a circle tangent to a straight 
string in the limit w — > 0. This familiar shape is indeed 
easily obtained with a nylon string, a metal string, or 
hair. An example with a Silica nanowire is presented in 
FiglU The relation between tension and excess length 
fi* = 2ir 2 C/£ 2 for the open trefoil knot was in fact al- 
ready used as an ansatz in Refi, where it was experimen- 
tally checked and used to evaluate the bending rigidity 
of actin. Here we show that it is an exact result in the 
limit w — > 0. Checking the n 2 dependence of fx* by vary- 
ing the knot type opens a novel and challenging line of 
investigations for experiments. 



B. Curvature energy of thick knots 

Two upper bounds for the energy of stiff knots can be 
derived from the limit of finite in. First, a general up- 
per bound directly follows from excluded volume effects. 
Indeed, we have seen in section lHll that k < 2/w. Hence, 

£ = — I dsn 2 < 2C*4r . (46) 
2 J w z 

which may be re-written in terms of the dimensionless 
normalized energy: 

Using l|15|) , one may also use Sid as an upper bound for 
the energy of a given knot. We do not know the precise 
value of Sid, which is dictated by the geometry of the 
ideal knot. Nevertheless, a lower bound for may be 
found from a recent conjecture based on results for lattice 
knots^i, which states that 

K id >bnl' 2 (48) 

where b is a positive constant. Using the Schwarz in- 
equality as in section ITV1 we find that: 

b 2 C 

Sid > -jn c - (49) 
which may also be written as > b 2 n c /2. 

C. Experimental shapes for stiff knots 



IX. DISCUSSION 

In the following, we shall make some remarks, and 
briefly mention some open issues related to the present 
work. 



A. On the limits w — > and £ — > oo 

Throughout the present paper, we have analyzed the 
limit w — > 0. Nevertheless, in a given experimental sit- 
uation, it is difficult to perform a variation of the width 
of the filament. Furthermore, we have assumed that the 
bending rigidity modulus does not vary with the width. 
But it often does, as e.g. in continuum elasticity^, where 
C ~ w~ A . A more natural limit for experiments would 
be to take C — > oo with fixed w. Both limits are equiv- 
alent. Indeed, the important point is that w/C — > 0, as 
seen e.g. in the inequality 



We have performed rudimentary experiments in order 
to corroborate the results of our Monte Carlo simulations. 
We do not look for quantitative measurements, but we 
rather aim for a qualitative confirmation of the theory 
and simulations. 

The experiments are performed with a plastic tube of 
width w = 1cm, and length C — 80cm. To close the tube, 
we have joined the two ends using a small stick inserted 
in both ends. This closure allows tangential matching as 
well as free local rotation of one end with respect to the 
other at the contact point. Therefore, the tube cannot 
store twist, and the torsion energy is neglected. These 
experiments are imperfect, and cannot be used for quan- 
titative purposes. For example, the small stick was less 
stiff than the tube, leading to a larger curvature at the 
junction. Moreover, the tube undergoes plastic deforma- 
tion, and it did not come back to a straight shape after 
the experiments. Furthermore, solid friction occurs at 
the contact of the tube with itself. Therefore, the curva- 
ture energy of the tube may not be fully relaxed. 

Despite these imperfections, the experiments qualita- 
tively confirmed the three types of shape already ob- 
tained, as shown on Fig^| Circular braids are obtained 
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FIG. 10: Photographs of closed knotted tube: 3i, 5i, T(5, 3) 
(torus knot), 4i, 62 in 4- fold cage configuration, 52 in 80 con- 
figuration. 



for knots with n = i: 3i and 5i, for which i = 2, and 
also the T(5, 3) torus knot with i = 3. Moreover, the 
expected 4-fold cage configuration is obtained for the 4i 
and the 5 2 . Finally, 80 configurations is found for the 52. 
These results are in perfect agreement with the results of 
the theory and the MC simulations of section |VTT1 



D. The restricted curvature problem 

Other minimization problems share similarities with 
the minimization of the curvature energy £. An exam- 
ple is the minimization of the total length of a knotted 
filament which exhibits a minimum possible radius of cur- 
vature R + . This constraint may result from an excluded 
volume effect which limits the bending angle between ad- 
jacent units in a polymer or a macroscopic chain. 

Let us assume that curvature k of a knotted filament 
cannot exceed 1/-R+, i.e. k < 1/-R+- Our aim here is to 
study the minimum possible length £t of this filament 
in the limit where w — > 0. First, from the inequality 
k = J dsn < C/R+ combined with J^J, we obtain a lower 
bound for the length of a knot with a restricted curvature: 



£ > 2imR + . 



(50) 



Once again, we want to avoid angular points and knot lo- 
calization, which exhibit infinite local curvature > 1/-R+, 
and braid localization is expected in the limit w — ► 0. 

We then use the same strategy as before to determine 
upper bounds. A first upper bound is found from the 
PBL configuration, defined in Fig|3] Because of the dis- 
continuous character of the constraint k < 1 / R+ , we can- 
not formulate the shape optimization problem with dif- 
ferential equations analogous to Eas. i|7l8[l to obtain the 
minimal shape and energy. We therefore use a simple 
ansatz for the geometry, where loops are formed by arcs 
of circles. The precise shape is defined in Fig^J The 
total length is: 



£ = 2n [i?i(7r + 26>) + 2R 2 t 



(51) 



where the free variables are R\, R 2 - The constraint of 
restricted curvature reads: 



i?i > R + ; i?2 > R+ ■ 
The angle 8 is given by the relation: 

R 1 cos(0) =-R 2 (l-cos(0)) 



(52) 



(53) 



with < 8 < 7t/2. The minimization of C with the 
constraints is straightforward^ and leads to R\ = R% = 
R + , so that 8 = 7r/3. Therefore, the minimum length of 
the PBL configuration within our ansatz is: 



7tt 

Cpbl = 2nR + — 



(54) 



A second upper bound for & is deduced from the cir- 
cular braid of radius R + , whose length is 2iriR + . Finally, 
we obtain a result similar to l|32|) : 



2imR+ < C) < 2ir min 
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R 



+ ■ 



(55) 



Similar conclusions are also drawn: Firstly, We conclude 
that the minimum length scales with the bridge number: 
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FIG. 11: Ansatz for the shape of the Point-Braid and Loop 
(PBL) configuration. The shape is made of arcs of circles of 
radius Ri, and i?2- The axis passing through the center of 
the circles goes through the contact point between two arcs, 
and defines the angle 9. 

£t ^ nR+. Secondly, the minimization problem is read- 
ily solved when n — i: we have £t = 2iriR + and the 
shape is a circular braid. 

The restricted curvature problem was recently ad- 
dressed numerically in Rcf. 27 for closed knots. As ex- 
pected from the above result, the circular braid configu- 
ration was found for some torus knots (with n = i = 2). 
A configuration similar to the 4-fold cage was also found 
for the 4i knot. Whether both problems always lead to 
similar geometries is still an open question. 

A second result of RefiSi is that a transition occurs for 
a finite value of the width w of the string. Does such 
a transition exist for stiff knots? Extensive simulations 
would be needed in order to answer this question. 



X. CONCLUSION 

In summary, we find that stiff knot mechanics crucially 
depends on the type of knot via a surprisingly simple 
quantity: the bridge number n. Indeed both the equi- 
librium energy of closed knots and the equilibrium ten- 
sion for open knots behave as n 2 . Up to now, the open 
trefoil knot is the only knot which has been studied in 
experiments. More complex knots have been studied in 
the fluctuation dominated regime for polymers such as 
DNA 6 . We hope that our results will motivate some nu- 
merical and experimental investigation of the mechanics 
of more complex stiff knots. 

As a second central result, we find that braid local- 
ization, which was checked here on simple knots, is a 
general and robust feature of the entanglements of stiff 
strings. We shall mention two possible consequences of 
braid localization. Firstly, the curvature energy -and 
thus braid localization- should be irrelevant for flexi- 
ble polymers. But it may be relevant for entanglements 
of semi-flexible polymers and fibersi&S. Imposing tan- 
gential contacts between the strings, braid localization 




FIG. 12: Excluded volume tubes of width wi and W2- 

questions the usual tube model for polymers, which is 
based on a lateral confinement due to simple crossings 7 . 
Secondly, braid localization also implies localization of 
friction and strain variations, and may have some impor- 
tant consequences on knot-induced polymer and filament 
break-up^. 

Our work opens a new line of investigation towards 
the understanding of the geometry and mechanics of stiff 
knots. But much yet remains to be done, and we would 
therefore like to conclude with a list of open questions: 

(1) In the present work, we have essentially analyzed the 
equilibrium state. But we intuitively expect the num- 
ber of metastable states to increase with the knot com- 
plexity. Can this be quantified? (ii) The question of 
the metastable states naturally leads to a second ques- 
tion: how can we study local stability (i.e. stability with 
respect to small perturbations)? (iii) We have studied 
the limit w — > 0. What happens for finite w? Here 
we have shown that the equilibrium energy £^ increases 
with w. The results of Buck and Rawdon 2 « on the simi- 
lar restricted-curvature problem suggest that qualitative 
transitions may occur at finite w. (iv) As shown in basic 
textbooks such as Refsd torsion usually plays an im- 
portant role in filament mechanics. How can we include 
it in the present theory? We hope to report along these 
lines in the near future. 

We wish to thank P. Peyla, Y. Colin de Verdiere, S. 
Baseilhac, C. Lescop, M. Eisermann, and H. Meyer for 
helpful discussions. YCdV pointed out the proof of ap- 
pendix Bl. 

APPENDIX A: MONOTONIC INCREASE OF 
THE EQUILIBRIUM ENERGY WITH w 

In this section, we show that the curvature energy of 
a knotted string with tubular excluded volume strictly 
decreases as the diameter w of the section of the tube 
decreases. 

To do so, let us consider a filament denoted as (1), of 
length C, and with an excluded volume tube of width w±. 
At equilibrium, the central line of the tube (1) is in the 
configuration r*(s). Let us denote the energy of this con- 
figuration as £*. We then consider another filament of 
length C, and of width W2, with wi <w\. If the filament 

(2) is in the configuration r^(s), the tube (2) is inside the 
tube (1). Therefore, there is no self-intersections or self- 
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contact for the tube (2). Thus, there is no interactions 
of the tube (2) with itself. Hence, following the results 
of section ITTT1 r*(s) is not the equilibrium configuration 
for the tube (2). Therefore, its energy £^ in the config- 
uration rj(s) is larger than the equilibrium energy 
of the filament (2). 

We have shown that £* 2 < £^ Ul when u> 2 < wi, for 
any W\ and tu 2 . As announced in the beginning of this 
section, £* decreases as w decreases -or equivalently 
d w £* > 0. 



APPENDIX B: ANGLE AND KNOT 
LOCALIZATION 



In this appendix, we show that the curvature energy 
of a curve in 3D diverges in the presence of angles, or 
knot localization. To do so, we will show the equivalent 
statement that, in the presence of an upper bound So f° r 
the energy: 



£ 



C 



dsn < So ■ 



(Bl) 



no angle or knot localization can be present. 

Let us consider a part of the curve r(s), running from 
si to s 2 . From the Schwarz inequality and IjBljl . one has: 



ds k < 



S2 \ 1/2 

ds 



< (as -si) 



1/2 



2So 
C 



S2 \ 1/2 

ds K 2 

l 

1/2 



(B2) 



Angles 



We observe that: 
|t(« 2 )-t(«i) 



ds d s t 



/ ds nn < dsn (B3) 

J si J si 



where (t, n, b) is the usual Frenet frame along the curve. 
Combing l|B2(l and {B3j, we obtain: 



|t(s 2 )-t( Sl )|<( S2 - Sl )^(-^ 



2Sr 



1/2 



(B4) 



Ea. (|B4(l shows that the tangent vector t is continuous 
along the curve. Therefore, no angle can be present. 



2. Knot localization 

If the part of the knot between si and s 2 is knotted, 
then there is a lower bound for the total curvature which 



is similar to Eq.©. But since the knot is an open one, 
one should use the modified bridge number defined in 
section rvTrfl which we denote as ni_> 2 . We thus have: 



ds K > 27T7ll^2 



(B5) 



Using this inequality together with l|B2|) . we obtain 



1 {2So 



1/2 



(B6) 



Knot localization means that there is a knot between S\ 
and S2 while (si — s 2 ) — > 0. But (|B6|) shows that, as 
(s 2 — Si) —> 0, one necessarily has ni^ 2 = 0. Since 
ni_> 2 = corresponds to an unknotted string, there is 
no knot localization. 



APPENDIX C: SOME RELATIONS FOR A 
STRUCTURE OF LINE-BRAIDS 



Closed knots 



a. Relation between £ and fi 



We here derive a relation between S, L, and [i. Inte- 
grating (|21|l over the pth line-braid, we obtain: 



S° - Li°m p £° p = A p .[r] p 



(CI) 



where [r) p is the difference between r at the beginning 
and at the end of the line-braid. Summing I|C1|) over p, 
we obtain: 



£°-/i £-^A p .[r] p = ^r B . ]T A p = (C2) 

p=l B p->B 



where the last equality follows from (|24|l . Finally, one 
has: 



M = 



£ 
C 



(C3) 



This equality is true for all structures which obey I|21I22|I , 
with the boundary conditions 124125(1 . 



b. Expression of e 

We now derive a relation between the total normalized 
energy, and the normalized energy of its parts. Eliminat- 
ing [i between Eas. (|Cl|) and i|C3() . and rewriting energies 
in terms of the normalized energies e, we find 



e;-A p .[r] 



p m p C 



1/2 



01/2 rn p Cp 



(C4) 
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where e p = £ p C p /m p C. Summing over p, we get: 

N , N 1/2 

TOpC 



(a) 



p=l 



(C5) 



There are three interesting cases in which the second term 
in the parenthesis vanishes. The first one is the situation 
where line-braids are arcs of circles, for which A p = 0. 
In the second case, all line-braids are loops (i.e. they 
start and end at the same point), implying [r] p — 0. The 
third case is the case where all line-braids have identical 
energies and length. If there is N line-braids, /i = £ /£ = 
N£ p /Nm p £p = Sp/rripCp. Such an equality, combined 
with 1|C1|) . implies that A p .[r] p = 0. In these 3 cases, one 
finally has: 



[E- P s 1/2 



(C6) 



and from Eq.jnU), one finds Ea. lP|) . 



2. Open knots 

a. Relation between £ and C e 

In the case of open knots, the relation (|C1(I is still valid 
for each line-braid. But the sum in (|C2|) does not vanish, 
and we obtain: 

N 

£°- t ,°£° =^A,[r] p = £r B . £ A p = [r.A]j(C7) 

p=l B p^B 

where the index ± indicate the two open ends at z± — > 
±oo. Using Ea. l|2(J|) . we obtain 




8=71/1 




e=o 



x 
-o- 



e=o 



FIG. 13: Schematics of (a) the loop solution. The loop part of 
solution used in the text starts and ends at the place indicated 
by the black dot. (b) the 8 solution. 



Following the same lines as in the previous paragraphs, 
one finds that: 

,C10) 

« = i(5V + A,IH.-S ! (en) 



m p C 



^p I f 

\ ' 

Summing these relations over all line-braids lead to 



1 + el^ 



N 



= J2 m p ( e p- a p-Wp 



r 



mr,C 



1/2 



Y&12) 



APPENDIX D: ENERGIES OF SPECIAL 
CONFIGURATIONS 

1. Energy of the planar loop 



[r.A]t -> -fi°z+ + n°z- = —fx / dz, 



so that (|C7|I can be re-written as: 



where C e is defined in Ea.l|lD|l. 



(C8) In this appendix, we determine the energy of one pla- 

nar loop in the configuration of Fig|3^. For planar so- 
lutions, A is parallel to the plane, and Eas. (|7l8|l reduce 
to: 



(C9) 



b. Expression of e e 

Let us denote with the index + and — the parts of the 
structure related to the open ends. Their excess lengths 
(/ ds — J dz along these parts) are respectively denoted 
as £ e + and £ e -. The + part extends from z + — > +oo to 
the first contact point, whose abscissa along z is denoted 
as zo+. For the — part, zo- is defined in a similar way. 
We then define Aq = zq+ — zq-- 



\{d s ef -^0* + a cose) = o. 



(Dl) 



where A = |A|, and 9 is the angle between the tangent 
vector t and A. 

The abscissa x, defined in Fig^J is written as: 



where d s 9 is given by Ea. (|Ul|l . Using the variable change 
v = — (x cos , with a = A/fi, we find: 



-x(C 
x = a — 



1/2 ,a 



v dv 



fij J v (a 2 - w 2 ) 1 /2( w + 1)1/2 



(D3) 
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o (a 2 - v 2 )V 2 (v + l) 1 / 2 



As seen from an inspection of Figll3fa. the constraint f dv 

which selects the loop solution with a tangent vector 
along y at the boundary is: 

x{6 = n/2) - x(9 = 0) = 2(x(9 = tt/2) - x(0 max )),(D4) Substituting this value in Eq.fDHJ), we find e 



(D8) 



The numerical solution of Eq.{D5|) leads to a w 2.158. 



loop 



18.19. 



where 6 max is such that v max = — 1. This constraint 
leads to an equation for a: 



v dv 



1 ( a 2 _ V 2)l/2( w + 1)1/2 



= 0. 



/„ (a 2 -« 2 ) 1 /2( tJ + l)i/2 
Moreover, since k = \d s 0\, the energy reads: 

£ioo P = y fds(d s 8) 2 , 



and the total length is: 



C^loop ' — / 



(D5) 



(D6) 



(D7) 



These integrals are both re- written with the new variable 
v, and finally the normalized energy reads: 



o &loop£loop 
^loop 



_! (a 2 - u 2 ) 1 ^ 



(v + 1) 1 / 2 



„ dW \ a 2 _ v 2)l/2 



_! (a 2 -i; 2 ) 1 / 2 ( u + 1)V 2 



2. Energy of the 8 

For the '8', we use the same method as in the previ- 
ous section. The constraint is now the periodicity of the 
curve, which imposes 



(D9) 



where max is the maximum angle reached along the 
curve. This leads to the following condition after the 
change of variable to v: 



v dv 



1 (a 2 -i; 2 )i/2( u + 1)1/2 



0. 



(D10) 



which is solved numerically, and leads to a — 1.53. More- 
over, 



c 



dv 



(v + l) 1 / 2 

(q,2 _ „2)l/2 

dv 



_ x (a 2 -v 2 )V 2 (v + 1)V2 
We find numerically that e§ ~ 55.93. 
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does not necessarily mean that we have reached the global 
minimum. 

As expected, the orientation does not affect the physics. 
The reader interested in this point is invited to check the 
invariance of the results with respect to the variable change 
s — > — s. 

Defining a = R 2 /(Ri + R2), one has C = 2nR 1 (tt + 2(1 - 



a) _1 arccos(a). For fixed a, dn^C > 0. Therefore, the min- 
imum of C is reached for Ri equal to its minimum possible 
value, i.e. Ri = R+. For fixed R\, d a £ > 0. Therefore, a 
must be equal to its minimum possible value compatible 
with Ri — R + , i.e. a = 1/2. Therefore, R2 = R+. 



